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BIDIRECTIONAL REFLECTANCE MEASUREMENTS FROM AN 
AIRCRAFT OVER NATURAL EARTH SURFACES 

INTRODUCTION 

As part of a program to determine the reflectance patterns of solar radiation 
from natural surfaces, measurements were made over a variety of surfaces with 
a scanning radiometer mounted on a CV-990 aircraft. Of particular interest was 
the reflectance pattern dependence on surface type and solar zenith angle. Since 
many present meteorological satellites utilize narrow beam scanning radiometers 
and are not capable of measuring reflected energy in all directions, reflectance 
patterns for natural earth surfaces are needed to interpret these measurements 
in terms of the albedo of the earth. 


INSTRUMENTATION 

During the spring of 1966 and 1967 reflectance measurements were made 
with a Medium Resolution Infrared Radiometer (MRIR)* flown on NASA’s Convair- 
990 high altitude research aircraft. This radiometer was a model of the MRIR 
instrument flown on the Nimbus meteorological satellites. 

The MRIR is composed of five Casagrainian telescopes which focus incoming 
radiation reflected from a scan mirror into five optical channels characterized 
by different filters. Thermistor bolometers are used as detectors in all its 
channels. 

A photograph of the MRIR is shown in Fig. 1. The scan mirror is mounted at 
an inclination of 45° to the axes of the telescopes on the left, seen reflected in 
the scan mirror. The scan mirror rotates so that the 50 milliradian instantaneous 
field of view scans through a 360° arc in a plane perpendicular to the axis of ro- 
tation in 7-1/2 seconds. The aircraft-borne instrument, shown inside an aero- 
dynamic fairing in Fig. 2, scans in the vertical plane containing the longitudinal 
axis of the aircraft. Scanning is fore-to-aft in the downward looking direction, 
as depicted in Fig. 3. 


*The term "Medium Resolution (rifra-red Radiometer" (MRIR) is actually a misnomer, since the 
satellite instrument contains a channel which responds in part to visible radiation (.2 p. to 4/i), 
and the airborne instrument contains two such channels (0.55/z to 0.85 \±, and 0.2fi to 4/x). How* 
ever, the term "MRIR" has been so widely adopted as to compel its use here. 


Three MRIR channels are sensitive to infrared radiation. The other two, 
sensitive to solar radiation in the intervals 0.55-0.85/4 and 0.2-4 /t, were used in 
this study. Their effective spectral response functions are shown in Fig. 4. 


DEFINITION OF TERMS 

The geometric relationships involved in a study of cloud reflectances are 
shown in Fig. 5. Energy from the sun impinges on the earth’s surface at a zenith 
angle £ 0 and is reflected at a zenith angle <9, with an azimuth angle 0 measured 
with respect to the principal plane of the sun from the forward scattering direc- 
tion. The amount of energy reflected is generally dependent on these three angles. 
For this study, the reflectance zenith and viewing nadir angles can be considered 
to be equal.* 

The term bidirectional reflectance is often used to characterize the reflec- 
tion of abeam of solar radiation from a plane surface (as illustrated in Fig. 5). 
The effective bidirectional reflectance p can be expressed mathematically by 


£ 0 ) ~ -- - -- - — (ster" 1 ) (1) 

H* cos £ 0 

where N (0 , 0 , £ 0 ) is the effective radiance measured by a narrow beam radiom- 
eter, i.e., the integral over all wavelengths of the product of the actual radiance 
at the radiometer and the appropriate spectral response function shown in Fig. 4. 

H* is the effective solar irradiance and is obtained by integrating the product 
of the spectral solar irradiance, and the appropriate spectral response func- 
tion 0^ shown in Fig. 4 over all wavelengths. 


H* = H^dX. (2) 

Vo 

(The spectral solar irradiance given by Johnston 2 was used here.) 


*A beam emerging from the surface at a zenith angle of 90° is observed at a nadir angle of 87° 
from an aircraft flying at an altitude of 40,000 feet. The difference between the two angles de- 
creases to 0° as the viewing nadir angle goes to 0. 
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The effective directional reflectance, R ( £ 0 ), may be obtained by integrating 
the effective bidirectional reflectance over the upward hemisphere, once the 
dependence of p on 6, </>, and £ 0 is known. The effective directional reflectance 
can be expressed mathematically by 


t 




p{Q y 4> } £ 0 ) cos 6 sin Q d$ dq6 


( 3 ) 


From Eq. (3) it follows that the directional reflectance of an isotropically 
reflecting surface is ifp. 

Albedo is defined as the ratio of the amount of electromagnetic radiation 
reflected by a body to the amount incident upon it. 3 Measured earth surface 
albedoes, however, are dependent on the time and spectral interval of the measure 
ments. For global heat balance studies, Raschke defines albedo as ”. , , that 
part of the daily incoming flux of solar radiation which is reflected back to space. 4 
The daily time-averaged effective albedo A may be found by integrating the 
measured effective directional reflectances with respect to time from sunrise 
(t j) to sunset (t 2 ), 

f' 2 - 

R(J 0 (t) cos £ 0 (t) dt 

A = _ (4) 

r* 2 

I cos £ 0 (t) dt 


The albedo of a constant surface may be obtained from narrow beam satellite 
measurements, N(0, <£, £ 0 ) once models for the spectral and time dependent 
functions are made. 


EXPERIMENTAL PROCEDURE AND DATA REDUCTION 

To measure effective radiance, for effective directional reflectance calcula 
tions with the aircraft-borne MRIR, it is advantageous to take measurements at 
as many nadir angles and azimuthal directions as possible. Since the solar 
zenith angle and cloud patterns are constantly changing, it is necessary to 
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measure the effective bidirectional reflectance of a surface in as short a time 
as possible. The CV-990, however, requires ten minutes just to turn around to 
change azimuthal direction. The rosette flight pattern, shown in Fig. 6, was 
chosen as a compromise between requirements for a large number of azimuthal 
directions on the one hand and a fairly short period of time required to fly the 
pattern on the other. With this one-and-a-half hour flight pattern, effective bi- 
directional reflectances with nadir angles (0 £ 0 < 90°) are measured in five 
azimuths). By assuming symmetry around the principal plane, (directions 1 and 
2 in Fig. 6), effective bidirectional reflectar. es are obtained for eight directions. 
These reflectance measurements are then integrated over a hemisphere to obtain 
effective directional reflectance. 

The rosette pattern was changed in 1967, to check the assumptions of sym- 
metry around the principal plane of the sun, to that shown in Fig. 7. 

Each leg of a rosette lasted four to eight minutes, yielding thirty-two to 
sixty- four MRIR scans. These data were averaged by a computer in 1° nadir 
angle steps for each leg to produce one set of p(0, 4 > , £ 0 ) vs. 0 values for each 
leg. p(0, cjt>, £ 0 ) was integrated numerically from 0 = 0° to 0 - 7t/2 for a given 
azimuthal direction by Eq. (3). Each integrated value was plotted against <£ and 
a curve was fitted through the points. The area under the curve was measured 
with a planimeter in order to find the effective directional reflectance. 


RESULTS 

Table I, (page 14) lists flight number, type of surface, aircraft and surface 
height, range of solar zenith angles over which a complete rosette was flown, 
and the calculated directional reflectance. 

Table II, (page 61), lists flight number, type of surface, aircraft and surface 
height, and the range of solar zenith angles over which measurements in the 
principal plane of the sun where made. 

The average bidirectional reflectances are plotted against nadir angle for 
each azimuth direction, and a polar plot for each complete rosette with isopleths 
of bidirectional reflectance is shown in Figs. 8-74. 

Clouds, (Figs. 8-37, and 57-65) show a strong forward scattering peak and 
a slight backs c altering peak with a minimum around the nadir. Measurements 
in the ,2-4/x band are generally lower than those in the .55-. 85 p region indicating 
that albedo measurements should not be taken in a narrow spectral interval but 
rather, should be taken in a wide interval encompassing as much of the solar 
spectrum as possible. 
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The ocean surfaces, shown in Figs. 46-53, and 66-74, show a slight back 
scattering per*k with stronger forward scattering, The reflectance pattern is 
influenced by sun glitter in the forward scattering direction. 

A forest, Figs. 38-45, shows a pattern similar to the strato cumulus pattern 
with strong forward and backward scattering with a minimum at low nadir angles. 
Bidirectional reflectances observed over sandy desert soils show very little 
anisotropy, as can be seen in Figs. 54-56, A slight back scattering peak has been 
observed in sand and grassland for very low sun angles by Salomonson 5 and 
Coulson, et al. 6 The measurements presented here do not show a forward 
scattering peak, presumably because the sun is not low enough in the sky to 
produce it. 

The reflectance of most of the earth surfaces measured shows strong 
anisotropy, especially in the principal plane of the sun. The degree of anisotropy 
is dependent on the type of surface and generally increases with increasing solar 
zenith angle. 
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Figure 1 -Medium Resolution Infrared Radiometer (MRIR). 
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C 0 = SOLAR ZENITH ANGLE 
9 = NADIR ANGLE OF REFLECTED ENERGY 

<(> = HORIZONTAL ANGULAR DEPARTURE OF 990'S FLIGHT PATH OUT 
OF THE PRINCIPAL PLANE OF THE SUN 

SZY - PRINCIPAL PLANE OF THE SUN. 

Figure 5-Dicgram showing nadir angles, azimuth and solar zenith angles. The sun’s 
energy enters with a zenith angle of £ 0 an< ^ ' s r ®F!ccted with an angle of 6 in the •p 
azimuthal direction. 
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Figure 6-1966 rosette flight pattern* 
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Figure 8-Bidirectional reflectances of strato 
cumulus tops at 4,500 ft., bases at 2,000 ft. 
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45 - STRATO CU 4,500 FT, (BOTTOM 2,000 FT.) OVER 
PACIFIC OCEAN (AIRCRAFT AT 35,300 FT.) 

MAY 4, 1967, FIT. 01 DAY 126 
40- AZIMUTH* 180,0 
ZENITH - 71 .4 
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Figure 9-Bidirectional reflectances of strato cumulus 
tops at 4,500 ft., bases at 2,000 ft. 
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Figure 10-Bidirectional reflectances of strato cumulus 
tops at 4,500 ft. bases at 2,000 ft. 
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Figure 11-Numerical integration of p over 6 vs. azimuthal direction for strato 
cumulus, tops at 4,500 ft., bases at 2.000 ft. over Pacific Ocean. 
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Figure 13-Bidirectional reflectances of uniform 
strato cumulus over Pacific Ocean. 
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Figure 14-Bidirectional reflectances of uniform 
strato cumulus over Pacific Ocean. 
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STRATO CU AT 6,300 FT. OVER FOREST (AIRCRAFT AT 41,000 FT.) 
MAY 11, 1966, FLT. 14 DAY 131 
AZIMUTH « 180.0 
ZENITH = 59.9 
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STRATO CU AT 6,300 FT. OVER FOREST (AIRCRAFT AT 41,000 FT.) 
MAY 11, 1966, FLT. 14 DAY 131 
45- AZIMUTH = 270.0 
ZENITH = 68 .5 


( .55- .85ja) 


(,2-4p) 


NADIR ANGLE ( 0 ) 

Figure 1 7—B idirectiona I reflectances of strato cumulus over forest. 
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Figure 19-Numerical integration of p over Q vs. azimuthal direction 
for strato cumulus over forest plotted. 





: * 

A 

0 1 L 1 1 1— i 1 1 J 

0 10 20 30 40 50 60 70 80 90 

NADIR ANGIE (0 

i-: 

. . £ 

f 

1 

1 

: A \ 

Figure 20-Bidirectional reflectances of strato 
cumulus at 2,300 ft. over Pacific Ocean. 
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STRATO CU OVER 

PACIFIC OCEAN (AIRCRAFT AT 40,000 FT.) 
JUNE 10, 1966, FLT. 34 DAY 161 
AZIMUTH- 180,0 
ZENITH -46.0 



NADIR ANGLE (0) 

STRATO CU OVER 

PACIFIC OCEAN (AIRCRAFT AT 40,000 FT.) 
JUNE 10, 1966, FLT. 34 DAY 162 
AZIMUTH » 90.0 
ZENITH =52.7 
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Figure 21-Bidirectional reflectances of stroto cumulus 
at 2,300 ft. over Pacific Ocean. 
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STRATO CU OVER 
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JUNE 10, 196 6, FIT. 34 DAY 162 
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STRATO CU OVER 

PACIFIC OCEAN (AIRCRAFT AT 40,000 FT.) 
JUNE 10, 1966, FLT. 34 DAY 162 
AZIMUTH = 225.0 
ZENITH = 51 .1 
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Figure 22-Bidirectional reflectances of strato cumulus 
at 2,300 ft. over Pacific Ocean. 
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Figure 23-Numericol integration of p over 6 vs, azimuthal direction for strato 
cumulus tops at 2,300 ft, over Pacific Ocean. 
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Figure 24-Bidirectional reflectances of strato 


cumulus over the Pacific Ocean. 
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STRATO CU OVER 
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Figure 25-Bidirectional reflectances of strato 
cumulus over the Pacific Ocean. 
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STRATO CU OVER 

PACIFIC OCEAN (AIRCRAFT AT 25,000 FT.) 
JUNE 15, 1966, FLIGHT 37 DAY 166 
AZIMUTH = 45.0 
ZENITH = 72,2 
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STRATO CU OVER 

PACIFIC OCEAN (AIRCRAFT AT 15,000 FT.) 
JUNE 15, 1966, FLIGHT 37 DAY 166 
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Figure 26-Bidirectional reflectances of strato 
cumulus over the Pacific Ocean. 
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Figure 29-Bidirectional reflectances of strato cumulus at 4,500 ft., 
bottom at 2,000 ft., over the Pacific Ocean. 
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STRATO CU 4,500 FT. (BOTTOM 2,000 FT .) OVER 
PACIFIC OCEAN (AIRCRAFT AT 25,000 FT.) 

MAY 6, 1967, FLT. 01 DAY 126 
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ZENITH = 65.7 
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Figure 30-Bidirectional reflectances of strato cumulus at 4,500 ft 
bottom at 2,000 ft., over the Pacific Ocean. 
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Figure 31 -Numerical integration of p over 0 vs. azimuthal direction for strato 
cumulus tops at 4,500 ft., bases at 2,000 ft. over the Pacific Ocean. 
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Figure 32-Qidirectional reflectances of broken 
strato cumulus at 15/000 ft. over Pacific Ocean. 
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Figure 33-Bidirectiona| reflectances of broken strato cumulus 
at 15,000 ft. over Pacific Ocean. 
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Figure 34-Numerical integration of p over 8 vs, azimuthal direction for 
broken strato cumulus tops at 10,000 ft, ovsr Pacific Ocean. 
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Figure 35-Bidirectional reflectances of strato 
cumulus at 10,000 ft. over Pacific Ocean. 
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Figure 36-Bidirectionoi reflectances of strato cumulus 
at 10,000 ft. over Pacific Ocean. 
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BROKEN STRATOCU AT 10-15,000 FT. OVER 
PACIFIC OCEAN AIRCRAFT AT 35,000 FT. (.55-. 85 p) 
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Figure 37— Numerical integration of p over 6 vs. azimuthal direction for broken 

strato cumulus tops at 10-15,000 ft. 
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PATCHES OF STRATUS OVER FARM & WOODED AREAS 
(AIRCRAFT AT 33,000 FT.) 

JUNE 3, 1967, FLIGHT 11 DAY 154 
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Figure 39-Bidirectional reflectances of farmland and 
wooded areas with patches of stratus. 
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Figure 40~~Bidirectional reflectances of farmland and 
wooded areas with patches of stratus. 
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PATCHES OF STRATUS OVER FARMLAND & WOODED AREAS 
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Figure 41 -Numerical integration of p over Q vs. azimuthal direction of 
farmland and wooded areas with patches of stratus. 
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DESOTO NATIONAL FOREST (AIRCRAFT AT 37,000 FT.) (,55-.85*j) 
SOLAR ZENITH 51 .6-68.6° MAY 10, 1966 
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Figure 45-Numericol integration of p over 6 vs. azimuthal direction 
for the Desoto Nationas Forest. 
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Figure 46-”Bidirectioncvl reflectances of the Pacific 
Ocean at ?*'/° 40 ' N 122° 48 W. 
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Figure 47— Bidirectional reflectances of the Pacific 
Ocean at 37° 40' N 122° 48' W. 
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Figure 48-B (directional reflectances of the Pacific 
Ocean at 37°40'N 122°4S'W. 
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Figure 49-Numerical integration of p over 6 vs. azimuthal direction 
for the Pacific Ocean at 37° 40' N 122°48'W. 
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Figure 51— Bidirectional reflectances of the Gulf of Mexico. 
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Figure 52-Bid irectiona I reflectances of the Gulf of Mexico. 
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TABLE H 


Fit. No. 

Type of Surface 

Surface 

Alt. 

Aircraft 

Alt. 

Range of 
Solar Zenith 
Angle (£ 0 ) 
[degrees] 



iO 3 

ft. 

1966: 13 

Strato Cumulus Over 
Okeefenokee 


40 

80.0 

85.6 

16 

Carson Sink 

3.9 

39 

40,7 

41.6 

32 

Salt Lake Desert 

4.3 

40 

21.2 

22.5 

34 

Strato Cumulus Over 
Pacific Ocean 

2.3 

2.5 

26.8 

28.8 

34 

Strato Cumulus Over 
Pacific Ocean 

2.3 

8 

30,5 

32.5 

34 

Strato Cumulus Over 
Pacific Ocean 

2,3 

16 

32.7 

33.7 

34 

Strato Cumulus Over 
Pacific Ocean 

2.3 

24 

35.8 

38.3 

34 

Strato Cumulus Over 
Pacific Ocean 

2.3 

32 

39.9 

40.8 

34 

Pacific Ocean 


16 

70.7 

72.6 

36 

Strato Cumulus Over 
Pacific Ocean 


40 

63.9 

73.3 

37 

Pacific Ocean 


40 

52.3 

50.4 

39 

White Sands 

4.0 

20 

9.6 

10.8 

41 

Pacific Ocean 


2.5 

25.5 

27.2 

41 

Pacific Ocean 


8 

22.2 

23.8 

41 

Pacific Ocean 


16 

18.5 

19.6 

41 

Pacific Ocean 


24 

15.5 

16.1 

41 

Pacific Ocean 


32 

15.1 

15.2 

41 

Pacific Ocean 


40 

22.0 

2.32 

1967: I 

Pacific Ocean 

4.5 

29 

85.9 

86.7 

I 

Pacific Ocean 


39 

88.2 

89.8 
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Figure 54-Bidirectional reflectances of carson sink, 
(aircraft at 39,000 ft) May 13, 1966, Day 133, Fit. 16. 
Solar Zenith 40.7-41.6°. 
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Figure 55-Bidirectional reflectances of Salt Lake Desert, 
(aii craft at 40,000 ft) June 9, 1966, Day 160, Fit, 32* 
Solar Zenith 22.5-21.2°. 
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Figure 56— Bidirectional reflectances of White Sand Desert, 
(aircraft at 20,000 ft), June 16, 1966, Day 167, Fit. 39. 

Solar Zenith 9.6-10.8°. 
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Figure 57-Bldirectional reflectances of Strato Cumulus 
at 1900 ft. over the Okeefenokee, (aircraft at 40,000 ft) 

May 10, 1966, Day 131, Fit. 13. Solar Zenith 80.8-85.6*. 
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Figure 58-B'idirectiona I reflectances of Strato Cumulus i 

at 1900 ft, over the Okeefenokee, (aircraft at 40,000 ft) I 

May 10, 1966, Day 131, Fit. 13. Solar Zenith 80.8-85.6°. I 

l 

I 




AVERAGE BIDIRECTIONAL REFLECTANCE ((>) AVERAGE BIDIRECTIONAL REFLECTANCE (p) 



Figure 59-Bidirecfionai reflectances of Strato Cumulus at 2,300 ft. 
over Pacific, (aircraft at 2,500 ft) June 10, 1966, Day 161, Fit, 34. 
Solar Zenith 26.8-28.8°. 
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Figure 60-Bidirectional reflectances of Strato Cumulus at 2,300 ft. 
over Pacific, (aircraft at 8,000 ft) June 10, 1966, Day 161, Fit. 34. 
Solar Zenith 30.5-32.5°. 
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Figure 62— Bidirectional reflectances of Strato Cumulus at 2,300 ft. 
over Pacific, (aircraft at 24,000 ft) June 10, 1966, Day 161, Fit, 34. 
Solar Zenith 35.8-38.3°. 
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Figure 63-Bidirectional reflectances of Strato Cumulus at 2,300 ft. 
over Pacific, (aircraft at 32,000 ft) June 10, 1966, Day 161, Fit. 34. 
Solar Zenith 39.9-40.8°. 
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Figure 65- Bidirectional reflectances of Strata Cumulus over Pacific 
Ocean (aircraft at 40,000 ft.), June 15, 1966, Day 166, Fit. 36. 
Solar Zenith 63.9-73.3°. 
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Figure 66-Bidirectional reflectances of Pacific Ocean at 39° N 
127° W (aircraft at 40,000 ft.), June 15, 1966, Day 166, Fit. 37. 
Solar Zenith 52.3-50.4°. 
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Figure 67-Bidirectional reflectances of Pacific Ocean at 38° 30 ' N 
126° 30 ' W (aircraft at 2,500 ft), June 17, 1966, Day 168, Fit. 41. 
Solar Zenith 25.5-27.2 
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Figure 68-Bidirectiondl reflectances of Pacific Ocean at 38° 30' N 
126° 30 'W (aircraft at 8,000 ft), June 17, 1966, Day 168, Fit, 41. 
Solar Zenith 22.2-23.8°. 
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Figure 69-Bidirectiona! reflectances of Pacific Ocean at 38° 30' N 
126° 30' W (aircraft at 16,000 ft), June 17, 1966, Day 168, Fit. 41. 
Solar Zenith 18.5-19.6°. 
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Figure 74— Bidirectional reflectances .of Pacific Ocean at 34 45 N 
121° 15' W (aircraft at 39,000 ft), May 6, 1967, Day 126, Fit. 1. 
Solar Zenith 88.2-89.8°. 




